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ABSTRACT 



r-| I We estimate the temporal change of magnetic flux perpendicular to the solar 

O-i' surface in a decaying active region by using a time series of the spatial distri- 

O ' bution of vector magnetic fields in the photosphere. The vector magnetic fields 

'^ ' are derived from full spectropolarimetric measurements with the Solar Optical 

"^ ■ Telescope aboard Hinode. We compare a magnetic flux loss rate to a flux trans- 

port rate in a decaying sunspot and its surrounding moat region. The amount 
>■ ■ of magnetic flux that decreases in the sunspot and moat region is very similar to 

^ ■ magnetic flux transported to the outer boundary of the moat region. The flux 

cn ' loss rates [{dF/dt)ioss] of magnetic elements with positive and negative polarities 

are balanced each other around the outer boundary of the moat region. These 
O ' results suggest that most of the magnetic flux in the sunspot is transported to 

the outer boundary of the moat region as moving magnetic features, and then 
removed from the photosphere by flux cancellation around the outer boundary 
of the moat region. 
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INTRODUCTION 



How and where is magnetic flux in sunspots removed from tlie pliotosphere? Sunspot 
umbrae are sometimes split by forni ation of a light bridge, which is a bright lane crossing 
the umbra (JBray &: Loughheadlll964j ). Smal l magnetic features with a typical size less than 
2" called moving magnetic features (MMFs; iHarvey fc Harveylll973[ ) are generally observed 
in the moat region that surrounds a decaying, mature sunspot. It has been reported that 
MMFs appear no t only in the decaying phase of sunspots but also in the growing phase (e.g. 
Wang et al.lll99ll ). The MMFs mostly appear around the outer boundary of the sunspot, 
moving almost radially outward during their lifetime ranging from a few minutes to 10 hr 
( jHarvey fc Harveyl 1 19731 : IZhang et al.l l2003l : iHagenaar fc Shind 120051 ) . The formation of the 
light bridges and MMFs is closely related to the fragmentation and disintegration of the 
sunspot magnetic flux. Indeed, it has been observed that the net fl ux carried away frorn . 



the sunspot by M MFs is larger than the flux decrease in the sunspot (JMartmez Pilletl 12002 



Kubo et al.l 120071 ) . This indicates that MMFs can be responsible for the flux loss of the 
sunspot. 

The mutual apparent loss of magnetic flux is often observed in the line-of-sight mag- 
netograms when one magnetic polarity element collides with another polarity magnetic 
element in the photosphere. This apparent flux loss is called "magnetic flux cancella- 
tion" as a descriptive term. It is observed that moving magnetic features often collide 
with apparently static opposite polarity magnetic features around the outer boundary of 



the moat region (IMartin. Livi. &: Wang Il985l : lYurchyshyn &: Wangl l200ll : IChae et al.l 12004 



Bellot Rubio fc Beckl |2005| ) and widely believed that understanding the flux cancellation 



process around the moat boundary is the key to understand the dissipat i on of the sunspot 
flux from the photosphere. Three models have been proposed by IZwaaru (119871 ) to describe 
the flux cancellation: (1) retraction of magnetic fields that connect an emerged bipole, (2) 
submergence of r2-loop formed by magnetic reconnection between the canceling two bipoles 
above the photosphere, and (3) emergence of U-loop due to reconnection below the photo- 
sphere. As expected in these models, h o rizontal magnetic fields are formed betwee n canceling 
magnetic elements (IWang fc Shilll993l : IChae et al.ll2004j : iKubo &: Shimizul 120071 ). However, 
whether upward or downward motions are obser ved in the cancellation sites depends on th e 
events and the positions in the cancellation sites (JHarvey et al.lll999l : lKubo fc Shimizull2007l ). 
Therefore, the nature of the physical process driving magnetic flux cancellation is still an 
open issue. 

This study attempts to address a basic question how much magnetic flux is carried away 
from the sunspot to the outer boundary of the moat region and is subsequently removed from 
the photosphere. Because it has been difficult to measure the magnetic field vector under 
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stable seeing conditions for the period longer than a typical lifetime of MMFs, the flux loss 
rate of the sunspot, flux transport rate due to MMFs, and flux cancellation rate have been 
independently estimated by using different data sets. A time series of spect ropolarimet- 



ric measurements with the Solar Op tical Telescope (SOT; iTsuneta et al.ll2008l ) aboard the 



Hinode satellite (JKosugi et al.l 120071 ) allows us, for the first time, to estimate an accurate 
flux change without any effects of atmospheric seeing. Moreover, the high spatial resolution 
observations with the SOT decrease the likelihood of spurious magnetic cancellation events, 
i.e., those for which magnetic elements with opposite polarities are located entirely within a 
resolution element and will dramatically increase the reliability of the results presented. 



2. OBSERVATIONS 

NOAA AR 10972 emerged on 2007 October 5 and formed a small bipolar sunspot, 
as shown in Figure [TJ Both the leading and following sunspots completely disappeared at 
the end of October 8, leaving only plage regions. The Hinode SOT started to observe this 
active region from 14:00 on October 6, so that we did not observe the growing phase of the 
active region. However, SOT provided a good data set for the decay phase. We selected 
observations of the following sunspot from 15:05 on October 6 to 07:12 on October 8, at 
which time the sunspot had significantly decayed (Fig. [2]). The reason why we focused on 
only the following sunspot was that a part of moat region surrounding the leading sunspot 
was outside the field of view. In this period, the spectropolarimeter (SP) of the SOT scanned 
the active region with the field of view of 152" x 164" every 1-2 hr, except for two 4 hr gaps. 
This provided us a time series of spatial distributions of the full polarization state for two 
photospheric Fe I lines at 6301.5 A and 6302.5 A. The slit scanning step was 0.297" with 
an integration time of 3.2 seconds (Fast mapping mode). SP took about 32 minutes to 
complete each scan. The pixel scale along the slit was 0.320", which was binned 2 pixels 
with the original spatial resolution of SP. In the same period, the narrowband filter imager 
(NFI) of the SOT continously obtained Stokes / and V images of the active region at a 
wavelength of -172 mA from the center of the lower chromospheric Na D line at 5896 A with 
a 2 minute cadence. The field of view was 276" x 164" with a pixel sampling of 0.16". 
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3. DATA ANALYSIS 

3.1. Magnetic Flux Density 

The temporal change of magnetic flux in and around the following sunspot was estimated 
from observations with the SP. The Stokes proflles of the two Fe lines were calibrated with 
a standard routine ("SP_PREP"; B. W. Lites et al. 2008, in preparation). The magnetic 
field vector and thermodynamic parameters were derived from the calibrated Stokes profiles 
with a Milne-Eddington Stokes inversion (T. Yokoyama et al. 2008, in preparation). The 
inversion code provided us magnetic field vector in a line-of-sight frame. A two-component 
model atmosphere, in which the photospheric atmosphere was composed of a magnetized 
atmosphere (polarized light) and a non- magnetized atmosphere (non-polarized light), was 
assumed in the inversion code. We estimated a fiux density of magnetic field vertical to the 
solar surface for the pixels that have the degree of polarization larger than 0.5 % as 

F = ^/|B|cos7, (1) 

cos 6* 

where a heliocentric angle {6), field strength (|B|), inclination angle (7) with respect to the 
local vertical, and filling factor (/). The filling factor is an areal percentage of the magnetized 
atmosphere in each pixel. 

A disambiguation of azimuth angles was needed to convert from the inclination with 
respect to the line-of-sight direction into the inclination with respect to the local vertical. To 
perform the disambiguation, we selected the azimuth angles closer to the azimuth of potential 
fields at the photospheric surface, and then interactively determined the azimuth to reduce 
discontinuities of azimuth arid incl ination angles by using the AZAM utility (written in IDL 



by P. Seagraves; iLites et al.lll995l ). The active region has a simple bipolar structure with 



simple unipolar spots. For such cases the ambiguity resolution with the AZAM is almost 



successfully performed (iMetcalf et al.ll2006l ). Furthermore, the sunspot was located near the 
disk center and the heliocentric angle ranges from 11° to 25°. This means that projection 
effects of errors in magnetic field inclination due to the disambiguation of azimuth angle are 
small. 

We applied an image cross-correlation for the magnetic fiux density maps in order to 
align the SP maps obtained at different times. The area around the following sunspot was 
used in the image cross-correlation in order to accurately remove the proper motion of the 
following sunspot. We did not use continuum intensity maps in the image cross-correlation 
because magnetic features were still observed after the following sunspot became very small 
and the granulation patterns had changed significantly. This allows us to trace the change 
of magnetic fiux until the disappearance of the sunspot. 
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3.2. Horizontal Velocity of Magnetic Elements 

Horizontal velocities of magnetic elements are necessary for estimation of a flux transport 
rate. Hereafter, the Na D line-of-sight magnetogram is deflned as the Stokes V image divided 
by the simultaneous Stokes / images. For each SP map (about 32 minutes), 16 line-of-sight 
magnetograms with a 2 minute cadence were obtained. We made 8 horizontal velocity maps 
with a 4 minute cadenc e from the 16 line-of-sight magnetograms, using a local corr elation 
tracking method (LCT; [November fc Simon! Eossl : Ichae et all I2OO1I : ISakamotol 12004 ) . The 
apodization window was a Gaussian with the FWHM of 1" in the LCT. The LCT was 
applied for the pixels that have Stokes V/I signals larger than 0.0015 in two sequential 
line-of-sight magnetograms. The threshold of 0.0015 was about la noise of the line-of-sight 
magnetograms. We assumed that the la noise of the line-of-sight magnetograms corresponds 
to the width parameter (the standard deviation) of a Gaussian fltted to a central part of 
a histogram for the signals in the line-of-sight magnetogram (Fig. [3]a). When a pixel had 
a coss-correlation coeflicient less than 0.9 or the Stokes V/I signal less than 0.0015, the 
horizontal velocity for the pixel was assumed to be km s~^. 



Finally, we averaged over the 8 horizontal velocity maps, and then aligned the averaged 
map to the SP map obtained at the same period, using the image cross-correlation between 
the line-of-sight magnetic flux density with the SP and the line-of-sight magnetogram with 
the NFL Figure[3?' shows a histogram of the horizontal velocity. The average of the horizontal 
velocity is 0.5 km s~^ , which is similar to the ave r aged speed of MM Fs and moat flow in 
previous studies (e.g. iBrickhouse fc Labontd llQSSi : IZhang et al.l 120031). Magnetic elements 
that move faster than the averaged moat flow of 0.5-1.0 km s~^ (JHagenaar fc Shind l2005l ) 



are also detected. Furthermore, radial outward motions of magnetic elements can be seen 
around the following sunspot in Figure [3]c. Thus, radial outward motions of MMFs are 
successfully obtained with the LCT. 



RESULTS 



Many moving magnetic features (MMFs) are observed outside the sunspot (Fig. ^). 
The MMFs of positive polarity (the polarity of the following sunspot) are dominant in the 
region with a radial distance (r^) less than 20" from the center of the following sunspot. 
MMFs of both polarities are located in the region r^ > 20". Most magnetic elements with 
negative polarity are located in the north-western side to the sunspot, and are in contact 
with positive polarity elements. Hereafter, we call regions with r^ = 0" — 7", r^ = 7" — 20", 
and Tg = 20" — 40" as the sunspot region, the unipolar region, and the mixed polarity region, 
respectively. In regular decaying sunspots, the unipolar region corresponds to the inner and 



-6- 



middle moat region, and the mixed polarity region corresponds to the region around the 
outer boundary of the moat region. We investigate the temporal change of magnetic flux in 
these three regions, and also investigate how much magnetic flux passes the outer boundaries 
of the three regions. 

Positive magnetic flux decreases at a constant rate in the sunspot region, as shown in 
Figure HI In the unipolar region, the positive magnetic flux increases from 01:00 to 06:00 
on October 7 as a result of the flux emergence indicated by the arrows in Figure [21 The 
emerging bipoles appear close to r^ = 20" and the opposite polarities associated with the 
emergence separate in the the regions r^ < 20" (r^ > 20") for positive (negative) flux. This 
flux emergence also causes increase of the negative magnetic flux in the mixed polarity region. 
After the flux emergence, the positive flux in the unipolar region and the negative flux in 
the mixed polarity region decrease at similar rates. On the other hand, the positive flux in 
the mixed polarity region increases during this period. The negative magnetic flux in the 
sunspot and unipolar regions is much smaller than the signed fluxes in the other regions. 
We estimate an increase and decrease rate of magnetic flux in each region by a linear fit to 
the time profiles from 06:15 to 21:26 on October 7 in Figure HI The magnetic flux changes 
nearly linearly during this period. The result is shown as "rfF/rft" in Figure O Even when 
the sunspot becomes too small on October 8, the positive flux in the sunspot region still 
decreases at a constant rate. On the other hand, the positive flux does not show any increases 
and decreases in the unipolar region and mixed polarity region on October 8. This tendency 
already may be seen in the later part on October 7. The change of the tendency suggests the 
possibility of different flux loss process in remnant active regions, but further, more detailed, 
investigations of this phenomena are needed. 

Figure [6] shows temporal evolution of total magnetic flux at each radius from the sunspot 
center. The areas that have large positive magnetic flux move away from the sunspot toward 
the mixed polarity region. Such outward motion of the positive flux is due to magnetic flux 
transported by the MMFs of positive polarity. The outward motion of the positive flux stops 
around 30" from the sunspot center. Figure [2] also shows that most of the positive magnetic 
elements do not go out of the mixed polarity region, except for the positive elements located 
in the north-eastern side to the sunspot. Most of the negative magnetic flux is located around 
30" from the sunspot center, and converging motion of the negative flux is observed around 
the 30" from the sunspot center (indicated by the arrows in Fig. [6]). These results suggest 
that the boundary of a supergranular cell, which corresponds to the outer boundary of the 
moat region, is located around 30" from the sunspot center. 

We estimate a radial transport rate {Fy{rs)) of magnetic flux at each radius (r^) from 
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the sunspot center using the following formula: 

where F{rs) is the magnetic flux density of vertical fields (Eq. ([1])) and Vr{rs) is the radial 
component of horizontal velocity at r^. We set ri and r2 as r^il", respectively. The behavior 
of F^{rs) with time in Figure [7] shows that the motions of magnetic flux seen in Figure [6] (e.g. 
outward motion of the positive flux, converging motion of the negative flux) are obtained 
by using Equation ([2]). The inward flux transport rate for positive polarity elements around 
the sunspot center is probably incorrect. This anomaly is mainly due to the shrinkage of the 
sunspot, yielding spurious inward Vr{rs)'- the LCT does not work well when the contrast of 
line-of-sight magnetograms is low. However, the flux transport within the sunspot region is 
not the object of this study. We focus on the flux transport rates at the outer edges of the 
sunspot region as well as the unipolar and mixed polarity regions during the period when 
the flux change rate is calculated from Figure H] (06:15 to 21:26 on October 7). The flux 
transport rates averaged over this period are summarized as "Fj," in Figure [51 The direction 
of the arrows shows whether magnetic flux is transported radially inward or outward. 



5. DISCUSSION AND CONCLUSIONS 

The averaged flux change and averaged flux transport in Figure [5] are described in 
units of Mx s~^, so that we can compare these values directly. The observed flux change 
[{dF/dt)obs] in each region and the observed flux transport [(-F!u)ofes] at its boundaries would 
have a relationship: 

dF dF dF 

{—TTjObs = {—7-) Emerge ~ [—77) Loss =t [F^jobs- (3) 

The increase of magnetic flux due to flux emergence [{dF/dt) Emerge] can be neglected in this 
case, because we select the period without any significant flux emergence. Thus, we can 
estimate an actual flux loss rate [{dF/dt)Loss] in each region from the observations. 

The total of flux decrease rates in the sunspot and unipolar regions {dF/dt = -3.2 - 4.8 
= -8.0x10^^ Mx s~^) is almost equal to the flux transport rate at the outer boundary of the 
unipolar region for the positive polarity {F^ = 7.4x10^^ Mx s~^). This means that most 
of magnetic flux that disappeared in the sunspot and unipolar regions is carried away to 
the mixed polarity region. Note that we do not trace each magnetic element, and all of the 
MMFs that separated from the sunspot may not reach at the outer boundary of the unipolar 
region. However, we make Figure [5] from the observations for about 12 hr, which is twice as 



long as the period that magnetic elements with the average horizontal speed of 0.5 km s~^ 
need to move through the unipolar region of a width of 14". The increase of the positive 
flux in the mixed polarity region supports the migration of positive flux into the mixed 
polarity region. Both the increase of the positive flux in the mixed polarity region and the 
flux transport for the positive polarity elements at the outer boundary of the mixed polarity 
region are smaller than the positive flux transported from the unipolar region. Therefore, 
the magnetic flux that is carried away from the sunspot (and moat region) mostly disappears 
in the mixed polarity region, especially near the outer boundary of the moat region. 

One issue is that the positive flux carried away from the sunspot region {Fy = 7.8x10^^ 
Mx s~^) is bigger than decrease of the positive flux in the sunspot region {dF/dt = -3.2x10^^ 
Mx s~^). T his tendency was also reported in the previous work with a lower (about 1") spatial 



resolution (JKubo et al.ll2007l ). As a result of no flux emergence in the sunspot region, the flux 
transport rate should be less than the flux decrease rate in the sunspot region. That is to say 
that the flux transport rate is overestimated at the outer boundary of the sunspot region. 
In the calculation of horizontal velocities, we use the apodization window with 1", which is 
lower than the spatial resolution of the magnetic field maps. Such a lower spatial resolution 
of the horizontal velocity maps probably causes the overestimation of the flux transport rate 
at the outer boundary of the sunspot region. Fuzzy, small magnetic elements with a short 



lifetime have been obse rved around the outer boundary of decaying sunspots (IZhang et al. 



20071 ; iKubo et al.ll2008l ). These fuzzy magnetic elements have higher outward motion and 
smaller magnetic flux than those of usual MMFs. In the estimation of flux transport at the 
outer boundary of the sunspot, magnetic flux is mostly represented by MMFs with small 
horizontal velocity and large magnetic flux, but its horizontal velocity is represented by the 
fuzzy magnetic elements. We believe that these fuzzy magnetic elements correspond to a 
fluctuation of field strength or a fluctuation of inclination of penumbral magnetic fields, 
and thus do not contribute to the flux loss of the sunspot. Further investigation using 
spectropolarimetric measurements with a higher cadence will be necessary to know the nature 
of such fuzzy magnetic elements and their impact on the presented calculations. 

The magnetic flux of negative polarity decreases in the mixed polarity region, although 
the negative flux converges from the inner and outer boundaries of the mixed polarity region. 
Considering that the negative flux moves into the mixed polarity region with the average rate 
of 1.6x10^^ Mx s~^, the actual flux loss rate [{dF/dt)Loss] in the mixed polarity region may 
be as large as 3.9x10^^ Mx s~^ from Equation ([3]). This flux loss rate of negative polarity is 
balanced by the actual flux loss rate of the positive polarity (3.9x10^^ Mx s~^), which is a 
difference between the flux transport rate (F„ = 7.4 - 0.8 = 6.6x10^^ Mx s^^) into the mixed 
polarity region and the flux increase rate [dF/dt = 2.7x10^^ Mx s~^) there. The flux loss 
rates with both polarities in the mixed polarity region are consistent with the cancellation 
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rates in active regions (IChae et al.ll200d . |2004J : iKubo &: Shimizul 120071 ). Furthermore, most 



of the magnetic elements with negative polarity are located in contact with the positive 
elements. These results suggest that magnetic flux cancellation at the outer boundary of the 
moat region is essential for the removal of the sunspot magnetic flux from the photosphere. 
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Fig. 1. — (a) Continuum intensity and {b) magnetic flux density of vertical field in NOAA 
AR 10972. The continuum intensity is normalized to the mean intensity of the quiet area 
outside the sunspots. These panels are made from observations with the Hinode SP from 
15:05 to 15:37 on 2007 October 6. The box is identical to the field of view of Fig. [2l The 
positions in the vertical and horizontal axes are given with respect to the center of the solar 
disk. 
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Fig. 2. — Time series of (a) continuum intensity and (b) magnetic flux density of vertical 
field for the following sunspot in NOAA AR 10972. The solid, dashed, and dash-dotted 
circles indicate the outer boundaries of regions called the sunspot region (7" from the sunspot 
center), the unipolar region (20" from the sunspot center), and the mixed polarity region (40" 
from the sunspot center), respectively. The vertical and horizontal axes show the positions 
with respect to the disk center in units of arcseconds. The white arrows indicate an emerging 
bipole. 
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Fig. 3. — (a) Histogram of the magnetic signal in a line-of-sight NFI magnetogram (solid 
line). The gray line shows a result of a Gaussian fitted to a central part of the histogram. 
(&) Histogram of horizontal velocity calculated by the local correlation tracking technique 
that is applied for line-of-sight magnetograms. The horizontal velocities are averaged over 
the period taking each SP map (about 32 minutes), (c) Spatial distribution of a radial 
component of the horizontal velocity averaged over the period from 10:59:32 to 11:27:31 on 
2007 October 7 in units of km s~^. Positive corresponds to a radial outward motion. The 
circles are same as those in Fig. [21 
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Fig. 4. — Time profiles of total magnetic flux in the sunspot region (diamond), unipolar 
region (triangle), and mixed polarity region (cross) as determined from the SP maps. The 
black symbols show magnetic flux of positive polarity, and the gray symbols show magnetic 
flux of negative polarity. The time profiles in the period between the two dotted lines are 
used for calculation of the flux change rate (dF/dt) in Fig. [H 
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Fig. 5. — Summary of observations for the flux change rate [dF/dt] and the flux transport 
rate (F^) in units of lO^^Mx s~^ 
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Fig. 6. — Total magnetic flux at each radius from the sunspot center vs. time for (a) positive 
polarity and (b) negative polarity . The total magnetic flux with positive (negative) polarity 
represents the sum total of magnetic flux for positive (negative) magnetic elements that have 
the same radius from the sunspot center. The hatched areas with oblique lines represent 
the periods without any SP observations. The dashed, dash-dotted, and dotted lines show 
the outer boundaries of the sunspot region (7" from the sunspot center), the unipolar region 
(20" from the sunspot center), and the mixed polarity region (40" from the sunspot center), 
respectively. 
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Fig. 7. — Same as Fig. [6], but for showing the flux transport rate vs. time. The flux transport 
rate (F„(rs)) is calculated as Eq. ([2]). The white indicates magnetic flux transported away 
from the sunspot center. The value averaged along the white part of each line is identical to 
the flux transport rate (F^) in Fig. [51 



